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Abstract 

We study the behaviour of lattice momentum-space gluon prop- 
agators for a pure SU{2) gauge theory at finite temperature. We 
find out that the magnetic mass is 0.26(7^(T)T; we have repeated 
the same calculations in three dimensions. 
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1 Introduction 



The subject of the finite temperature behaviour of the gluon propagators 
has attracted attention due to its relevance in understanding processes in the 
early universe. Very early approaches |]T], have identified the two different 
kinds of masses arising when considering longitudinal or transverse degrees 
of freedom, called electric and magnetic masses respectively. It has been 
argued 0] that they are gauge invariant, thus measurable quantities. 

The pole for the electric gluon propagator may be obtained in perturba- 
tion theory, with the well-known one-loop result: 

me = ^ 6 3{T)T, (1) 

for Nf massless quark fiavours and the SU{Nc) gauge group. This has been 
measured numerically using the heavy quark potential and the correlations 
of Polyakov lines [|, |]. 

The magnetic gluon propagator has no (non-zero) pole at the one-loop 
level, however at two loops there is a contribution to the magnetic mass pro- 
portional to g'^{T)T. Higher loops contribute multiples of the same quantity, 
leaving no possibility for perturbative calculations. A non-zero value for this 
non-perturbative quantity may soften the severe infrared problems of finite 
temperature perturbative calculations, since it will act as an infrared cut-off. 
The subject of the magnetic mass has been attacked since the eighties, using 
various methods. On the lattice, it has been measured using the effects of 
twisted boundary conditions on bulk quantities |^ |^ and yielded the value 
0.2Ag'^{T)T. The next approach |^ has been more straightforward, in the 
sense that it studied the correlators of (gauge- variant) gluon operators. It 
found effective masses increasing with distance, being thus in conflict with 
the Kallen-Lehmann representation. The interpretation was that this odd be- 
haviour can be acceptable for a conflned state (such as the gluon), although 
it should be rejected for a physical particle. 

Restricting ourselves to pure gauge theories, we have three scales that 
enter the game: the temperature T, the scale of the electric mass, g{T)T, and 
the scale of the magnetic mass, g'^{T)T. At sufficiently high temperatures 
one expects that g{T) will be very small, therefore the electric mass will 
be much smaller than the temperature scale and the magnetic mass even 
smaller. This is not the case however in realistic lattice simulations, since 
at high temperatures the finite-size-effects become big, thus restricting the 
investigation to relatively low temperatures. 

There have been two recent determinations of the electric and magnetic 
masses; the first one measured correlators of gauge invariant objects in 
QCD, yielding rumag ~ (2.9 ± 0.2) T, ruei ^ (1.4 ± 0.2) T. The second one 
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1T0[] is based on the gluon propagators themselves in SU{2) gauge theory and 
yields rrimag = 0.466(15) g^{T) T, ruei = 2.484(52) T. 



On the other hand, there has been a different interpretation ||TT[ of the 
results of 0, claiming that the gluon propagators do not really behave as 
decreasing exponentials, but have instead a power law behaviour at small 
distances, which would mean that the results of are compatible even with 
a vanishing value for the magnetic mass. 

An approach that could shed some light on the space-time dependence of 
the gluon propagator is the study of momentum space propagators ([|12|, p!3|). 
This has the advantage that one studies the whole momentum space propaga- 
tor, rather than the few numbers which survive the sum over the hyperplanes. 
Moreover it opens the possibility to study the behaviour in each momentum 
region separately and relate the momentum under study to the scales that 
enter the problem, which is not possible in configuration space studies. Fi- 
nally, as shown in []T2[, there is the very important technical advantage of a 



much better behaved covariance matrix. In particular, the covariance matrix 
for configuration space propagators is singular for the lattice sizes presently 
used; this is not the case in momentum space, offering the possibility to 
perform fully correlated chi-squared fits. 



2 The finite temperature gluon propagator 

At finite temperature any second rank tensor may be expanded in the 
basis of the following four tensors: 

Pj. = - (2) 



Pi = (<S^-^)p(*--^) (3) 



pG ^ u\ 

p'^ = - ^) + (^^(^^^ - ^)) (5) 

In the Landau gauge the free gluon propagator is of the form: D^y{k) = 
fc^+k^ (-P^t, + -P/JL)- If one expands the self energy tensor II^i, on the above 
basis: 

uf^ih, k) = nfp^^, + u^'Pl + ufp'^^ + u^^p^,. (6) 

Taking into account that n}''(k, ki) = 6"'^Ilj(k, k^), where the index J may 
be L, T, S or G, one may show that the full gluon propagator at finite 
temperature is given by the expression: 
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Now we observe that Pj^{k) = -P^(fc) = 0, /i = 1, ...4 for any value of k, 
while Plj{k) = kl ^2% , hj = 1,2,3, which vanishes for all i, j when k^ = 0. 
This fact, that the static spatial self-energy is always transverse, may be 
shown |l[] using the Ward identities satisfied by making no appeal to a 
specific gauge. In the sequel we stick to the static case, k^ = 0. The limits 
of ni(k, A;4 = 0) and YlT^kyk^ = 0) as k goes to zero are the electric and 
magnetic masses (squared) respectively. In this work we are going to measure 
the (lattice versions of the) quantities 

C,(k) = lj2 6abGf^{k, k, = 0) = G{k, k, = 0)Pj:{k, k, = 0), (8) 

^ a,b=l 

where j runs from 1 to 3 and no sum over j is implied. In fact, we do not fit 
Cj(k), but the quantities 

The infrared non-perturbative behaviour of the gluon propagator has been 
studied using a variety of methods, with conflicting conclusions. Other works 
[llB] predict for the propagator a very singular behaviour at zero k, like 
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(k2) 

and a confining property is conjectured; on the other side, there are claims 
[|l^ that the tree-level pole at k^ = is removed, due to non-perturbative 
effects and the Green's function vanishes at zero momentum. It is a fact 
that in the infrared region there should appear a dynamically generated mass 
M{g, /i). However, the generation of this mass cannot be seen in perturbation 

constant 

theory, since it must behave like M{g,jj) ^ /ie for very small g, 

that is it should exhibit an essential singularity at zero gauge coupling. It 
is expected that non-perturbative effects may yield negative powers of the 
momentum in the vacuum polarization function of the gluon, introducing 
various mass scales: 

(10) 

g is the coupling constant and m{g, /i, T), b{g, /i, T) have dimensions of mass 
and depend non-analytically on g. We have only shown the first two terms of 
an expansion in j^, since they correspond to well-known suggestions about 
the gluon propagators. In particular, if b{g, /i, T) = 0, a non-zero m{g, /i, T) 
gives rise to a mass pole in the gluon propagator. On the other hand, if 
m{g, /i, T) = 0, a non-zero b{g, fi, T) will give rise to a propagator of the form 



which has been proposed by Gribov |[T5|. 
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3 The method 



We measured correlations of the gauge potential, defined through the 
relation: 

A^{x) ^ ' ^ = 1, 2, 3, 4. (12) 

The various expectation values should be calculated in a lattice version of 
the Landau gauge, implemented by making the quantity 

S[A] = Tr[A{n)U^{n)K\n + /x)] (13) 

a global maximum with respect to the gauge transformations A(n). This 
guarantees not only the satisfaction of the lattice version of the condition 
dfj^A^^ = 0, but also the additional constraint imposed by Gribov. Let us 
note that the lattice Faddeev- Popov operator is just the Hessian matrix of 
S[A] with respect to A[n]. The maximization condition fixes the sign of the 
Hessian to be the correct one. Moreover, the maximization condition also 
enforces the smoothness of the continuum limit. Of course, one cannot say 
whether this proposal really gets rid of all of the Gribov copies; however, it is 
easily seen that the Gribov copies of the trivial configuration are eliminated, 
since the latter is the only constant configuration which maximizes S[A]. 
Thus, it may be expected that also the problem in its general form may be 
less severe. 

In order to calculate the gauge dependent correlators, we updated with 
the usual Wilson action and transformed the resulting configurations to the 
Landau gauge before taking measurements. This can be proved [§] to in- 
corporate the effects of the Faddeev-Popov determinant. For SU(2), which 
we are considering, we sweep through the lattice and at each lattice site we 
calculate analytically the gauge transformation that maximizes the sum of 
the links beginning or ending at this site. Of course, this gauge transforma- 
tion will disturb the gauge condition on the neighboring sites, so we expect 
that the algorithm relaxes to the global maximum we are looking for after 
several sweeps. A good gauge fixing is very important for the reliability of 
the results. A straightforward check is to calculate the < As{x)A3{0) > cor- 
relator in configuration space, summed over the directions 1, 2 and 4: this 
quantity should be constant, as a consequence of the gauge fixing. We have 
considered the gauge fixing as good enough if the variation of this correlator 
was not larger than 0.1 percent. This in turn dictated that the quantity 
W]^^k,n,tM{d^A'^{n)Y (which must vanish in the Landau gauge) should be 
less than about 10~^. This last requirement has been practically used as the 
criterion to stop the gauge fixing iterations. Let us note that the number 
of measurements done varied between 3000 and 5000 for each point; we also 
mention that a number of overrelaxation sweeps ( 5) was performed between 
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successive Monte Carlo steps. In figure 1 we show the configuration space 
correlators < Ai{x)Ai{0) > and < A3{x)A3{0) > versus x. The first correla- 
tor decays as usual, while the second one is constant, because of the gauge 
fixing, as just explained. 



4 Results and conclusions 

Let us first explain the way of analyzing our results. We have measured 
the correlators of link variables Aj (j — 1,2,3) on the lattice and then 
performed the Fourier transform, considering only the case = 0. The 
lattice versions of the propagators have been considered in the fit function 
G(k^) = (yU = ma), which means that we used everywhere the di- 

mensionless lattice version k of the three-momentum rather than the usual 
one k. Let us recall the relevant definitions: pjO = ki = 2sin^ where 
kitt — ^ni,ni = — ^ + 1, . . . , Y for the even N we have been using. We 
have considered lattices with two different temporal extents {Nt = 2 and 
Nt = 3) to get some feeling about the finite size effects due to the smallest 
of the lattice dimensions. Of course the corresponding momentum k^ has 
been set to zero, as an external momentum, however it also appears as an 
internal momentum in the self-energy graphs and may very well influence the 
results. In most simulations the spatial dimension of the lattice has been 
10, however we have also used = 12 and = 14 in some cases. 

We have chosen to measure the momentum p {p = \p\, k = |k|) in units of 
the temperature T of each lattice, since this is one of the dominant scales of 
the problem (along with g{T)T and g^{T)T). We note that ^ = Ntkmax, 
so the range of the quantity explored by each lattice is proportional to 
its "temporal" extent, provided N is kept constant. Taking this fact into 
account, the (lattice) momentum range [0, kmax] has been divided into Nt 
)arts and the gluon propagators have been considered over the intervals: 
kmax^, kmax§-^ , n = 1, . . . , ATj, for the quantity k, corresponding to the 

intervals (n — i)kmax,'>T'kmax ^"^ — l)---)-^t for the quantity ^. We note 
that this division of the intervals is arbitrary; one might divide the intervals 
in a different number of subintcrvals (bigger or smaller than Nt). We observe 
that a specific value of n yields the same values of |; no matter which lattice 
we are considering. In particular, the value n=l corresponds to the results 
we are going to depict in the figures and has to do with the interval < 
y < kmax- Note that for the lattices 10^ x 2, at most ten momenta may 
be considered in the fits, while for the lattices 10^ x 3 the corresponding 
number of momenta is four. In particular for the Nt = 2 case we observed 
that the qualitative results we refer to in the following do not change if we 
consider shorter intervals within the ones corresponding to the various values 
of n; there are minor changes in the values of the mass and one gets smaller 
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values for the chis squared. There is a serious exception, though, for n = 1: 
the point k = should be included in the fit in any case, since the second 
smallest momentum possible {jj^) is already of the order of T. 

We have measured the propagators for several values of ranging from 
0.60 to 6.0 (the values greater than 4.0 could only be reached on the lattices 
with Nt = 2; for Nt = 3 the required Pg would be too large, given the spatial 
dimensions of the lattices we have been using). 

In figures 2a, 2b, we show the results for the spacelike propagators (mag- 
netic sector) as a function of the temperature for Nt = 2, 3 respectively, at 
"low" momenta (n=l). We plot the output for ^ = versus ^, as well 
as curves fitting these data . We have also included in the figures the data 
for ^ < 1. These have nothing to do with the magnetic mass, of course, 
however it is worth noting that no dramatic effects occur in the values of 
these screening masses at the phase transition point. One may observe that 
lattices with Nt = 2 and Nt = 3 give values quite close for equal ^ (with 
some small discrepancies around ^ = 2). 

To construct the fitting functions we invoke the expected behaviour m = 
cg'^(T)T of the magnetic mass along with the ansatz g~'^{T) = i^^og{^), 
inspired by (the leading order in) the renormalization group equation for the 
gauge coupling. In practice we have fitted — to the form A logiSr) + B and 

determined c, A through c = ^^^2^ , A = e~A. The fit has taken into account 
only the data with > 1.6. It turns out that 

c{Nt = 2) = 0.20(3), \{Nt = 2) = 0.55(15) 
clNt = 3) = 0.26(2), XiNt = 3) = 0.26(8). 

These values for c are more or less in agreement to the early publications 
on the subject, where c has been found 0.24. 

Let us note that using the above expression for g(T) and the just men- 
tioned values of A, we may find the ratio of the magnetic mass to the tem- 
perature: 

m{Nt = 2) = 1.9T, m{Nt = 3) = 1.6T 

For even larger momenta (n = 2, n = 3) the mass is zero in all cases. 
Thus, at high momenta the correlators tend to free field ones. 

It is of interest to check what happens to the quantities measured above in 
the case we have a three-dimensional lattice, which would correspond to the 
"infinite temperature" limit. In figure 3 we show the results for the resulting 
masses for a three - dimensional (12^) lattice. On dimensional grounds one 
expects that any mass should be proportional to g^, the square of the three- 
dimensional gauge coupling. We find that 

m = 0.49(3)^3^ 
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If we make the (crude) identification of the (bare) gl to g'^{T)T, we see that 
the three-dimensional result would imply a much larger value for the quantity 
c than the one found from the four-dimensional lattices; this result is, in fact, 



close to the result given in |]TU 



Let us note that, if we also treat 2x2 plaquette correlators in the same 
way as the link correlators, we get rupi = 2.3{5)g^ (in agreement with ||14||), 
that is a much larger factor than the one extracted from the link variables. 
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Figure captions 



Figure 1: The correlators < Ai{x)Ai{0) > (diamonds) and < A3{x)A3{0) > 
(crosses) versus a; for ^ = 2.5, N — 10 and Nt — 2. 

Figure 2: — versus ^ for the space-hke hnk correlators with n=l in the 
cases Nt — 2 (figure 2a) and Nt — 3 (figure 2b). The fitting curves are also 
plotted. 

Figure 3: ma versus G = g^a for the link correlators with n=l in the 
three-dimensional case. The fitting curve is also plotted. 
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